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In the previous paper, it was shown that holes doped into 
BaBiOs self-trap as small polarons and bipolarons. These 
point defects are energetically favorable partly because they 
undo locally the strain in the charge-density-wave (Peierls in- 
■ sulator) ground state. In this paper the neutral excitations 
of the same model are discussed. The lowest electronic exci- 
tation is predicted to be a self-trapped exciton, consisting of 
an electron and a hole located on adjacent Bi atoms. This 
excitation has been seen experimentally (but not identified 
as such) via the Urbach tail in optical absorption, and the 
multi-phonon spectrum of the "breathing mode" seen in Ra- 
man scattering. These two phenomena occur because of the 
Franck-Condon effect associated with oxygen displacement in 
the excited state. 
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I. INTRODUCTION 

The previous paper on BaBiOa [1| presents a simple 
and yet sufficiently realistic model |2| which gives a mi- 
croscopic description of the Peierls distortion, and makes 
unambiguous predictions about the low-energy behavior. 
The previous paper looks carefully at the small polarons 
created by doping or photoemission. These were noticed 
already by Yu, Chen, and Su Here we examine what 
happens when an electron-hole pair is inserted, and dis- 
cover a self-trapped exciton which has not previously 
been discussed. This excitation gives a natural expla- 
nation for the Urbach tail seen in optical absorption ^] 
and for the multi-phonon replicas of the oxygen breath- 
ing phonon (the Peierls mode) seen by Kozlov et al. ^] 
in infrared absorption, and by Tajima et al. |^ in Raman 
scattering. 

The model |^ , based on work by Rice and Sneddon , 
has nearest-neighbor hopping (with strength t) of elec- 
trons between Bi s-orbitals on a cubic lattice. Oxygens 
which lie between pairs of Bi atoms, experience a linear 
repulsive force g{n2 — ni) away from the Bi atom with 
the larger electron charge (ni) and toward the Bi atom 
with the smaller electron charge (7^2) • The oxygen atom 
also has a harmonic Einstein spring constant K. The 
model has a single dimensionless constant P = g^/Kt 
in adiabatic approximation (motion of the oxygen atoms 
neglected.) We choose parameters i=0.35 eV, K = 19 
eV/A^, g = 1.39 eV/1, yielding P = 0.30. 



II. ATOMIC LIMIT T = 

As in the previous paper, the zero bandwidth {t — 0, 
P — > 00) case provides a soluble limit which gives good 
insight into the more complete theory. If a hole is created 
in the lower Peierls or occupied valence band, and an 
electron in the upper Peierls or empty conduction band, 
without allowing any lattice relaxation, the energy cost is 
2A = 2Ag^ /K . The hole is a Bi'*+ ion on a site normally 
occupied by Bi'^"'", and the electron is a Bi^+ ion on a site 
normally occupied by Bi^+. 

Because of the simple Hamiltonian, the problem has 
particle-hole symmetry. This is true at all values of pa- 
rameters, not just i = 0. Upper and lower Peierls bands 



have energies ±A(A:) where \(k) — ^ IS? + e(fc)^, and e(fc) 
is the nearest neighbor "cosine" band dispersion. Thus 
particle and hole bands are symmetrical. 

If the electron and hole are spatially separated, each 
will lower its energy by lattice relaxation. Both electron 
and hole save 3 A /4 in elastic energy by moving surround- 
ing oxygens half way back to their cubic perovskitc posi- 
tions. However, each unpaired electron on a Bi^+ ion has 
its on-site energy increased by A/2, for a net lowering of 
energy by A/4 for each, or a cost of 3 A/2 to create the 
isolated pair of polarons. 

As shown in Fig. |l|, there is an advantage to having 
the electron and hole polarons form on adjacent Bi ions. 
The oxygen between the two Bi^+ ions can further re- 
lax to the undistorted perovskite position. Rather than 
12, now only 10 oxygens are half displaced, saving an 
additional A/ 12 = g^ / K in elastic energy. This bound 
pair of polarons, called a "self-trapped exciton," costs 
-Eoxc = 17A/12 rather than 2A to create. The total en- 
ergy gain from lattice relaxation is the exciton binding 
energy Ecxc = ^g^/if. Taking the on-site Coulomb in- 
teraction into account, the break-up of the Bi^+ paired 
electrons into separated electron-hole pairs will reduce 
the excited state energy further, to Ccxc — ^9^ /K + U 
where U is the Hubbard on-site Coulomb repulsion. 
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FIG. 1. Schematic structure of the self-trapped exciton. 
The electron is in the left-central cell ( "exciton-electron" ) , the 
hole ( "exciton-hole" ) in the right-central cell. Filled circles 
are Bi^^ ions with two s electrons, shaded circles are Bi*^ ions 
with one s electron, and open circles are Bi^^ ions with no 
s electrons. X's denote oxygen ions. The central oxygen ion 
is undisplaced, which lowers the energy relative to separated 
electron and hole polarons. 

The energy of the ground and excited state depend on 
the amount of lattice distortion. The configuration coor- 
dinate a smoothly connects the optimal oxygen configu- 
rations of the ground state at a = and the self-trapped 
exciton state at a = 1, as shown in Fig. |^. The en- 
ergy depends quadratically on a, separating ground and 
excited state by the gap 2 A at a = 0. 
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FIG. 2. Excited state and ground state energy versus con- 
figuration coordinate a for t — 0. At a = 0, the ground and 
excited state are separated by the gap 2A. The self-trapped 
exciton at a = 1 has its energy lowered by the maximal 
amount - the exciton binding energy toxc. Including on-site 
Coulomb repulsion, the energy curve of the excited state will 
be shifted down by the Hubbard U. 



III. NUMERICAL RESULTS 

For t we are not able to solve analytically, and rely 
on numerical studies. We find that the excited state has 
zero slope dE/da = at a = 0, and downward curvature. 
No potential barrier separates the excited state at a = 
and the self-trapped exciton at a = 1. This is contrary to 
most self-trapped exciton theories , where a metastable 
free excited state and a self-trapped exciton state can 
coexist, separated by a potential barrier. 

For t ^ Q, localizing and delocalizing terms in the 
Hamiltonian compete against each other, yielding a crit- 
ical coupling constant rc(exc) above which self-trapped 
exciton states exist. To obtain the lowest excited state 
numerically, oxygen positions are varied using the algo- 
rithm described in the preceding paper, under the con- 
straint that an electron was removed from the highest- 
lying valence band state and inserted in the lowest-lying 
conduction band state. We find rc(exc) « 0.175, close 
to the polaron stability limit Tc{P) = 0.18 of the previ- 
ous paper. Since the self-trapped exciton is an electron- 
polaron weakly bound to a hole-polaron, it is not sur- 
prising that the polaron and the self-trapped exciton 
have similar threshold coupling constants. The first cor- 
rection to strong coupling for eexc can be approximated 
using a "vacuum correction" of order 1/F^ discussed in 
Sect. in. A of the preceding paper: eexc/^ ~ 7F— l/2rand 
approaches asymptotically the i = limit eoxc = . 
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FIG. 3. Exciton binding energy toxc for t 7^ computed 
numerically. The bold solid line is the i = prediction. The 
dashed line is an approximate l/F'^ correction. 

The behavior of self-trapped exciton states and their 
evolution with F is very similar that of polarons |l[ . The 
local change in the oxygen environment leads to the ap- 
pearance of localized gap states (Fig.||) in addition to the 
states which constitute the exciton. 
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FIG. 4. Energies of localized single particle states which 
accompany the self-trapped exciton states lying in the gap. 
The spectrum has particle-hole symmetry. Five doubly oc- 
cupied states split from the valence band, corresponding to 
states surrounding the exciton-electron. The mid-gap states 
derived from the valence and conduction band respectively are 
occupied singly each, corresponding to the exciton-hole and 
exciton-electron. Five empty states split from the conduction 
band. The ten states near the conduction and valence band 
edges are localized on the 10 Bi atoms which surround the 
two Bi atoms which constitute the exciton. 

Six states emerge from the valence band and six from 
the conduction band, symmetrically around E = 0, be- 
cause of electron-hole symmetry. The lowest five states 
have a small upward shift {g'^/K when t = 0) from 
the top of the valence band, and remain doubly occu- 
pied. These are electrons on the five A sites surrounding 
the B sublattice site now occupied by the excited elec- 
tron (called the "exciton-electron"). Their degeneracy is 
partly lifted when they couple to other sites by hopping, 
shifting their energies toward the valence band. The next 
state up also derives from the valence band, and is oc- 
cupied once. It is the electron that remains behind on 
an A site when one electron is removed. Alternately, 
this is the site of the "exciton-hole". Its energy {Tg^/K 
above the valence band when t—Q) is slightly higher than 
for an isolated hole state {Qg^ /K) because one of the 
six surrounding oxygens, shared by the exciton-electron 
and the exciton-hole, is fully relaxed to the symmetric 
position. With decreasing coupling strength, the states 
spread out, as can be seen in the continuous evolution 
of their inverse participation ratios (IPR) in the inset in 
Fig. |5| . The IPR I /Pi for a state i is defined as 

l=^|vl/,(n)|4 (1) 

where the sum runs over all sites. It measures the lo- 
calization of a state. An IPR of 1/P represents a state 
that is localized on P atoms. The mean radius (r) of 
the two states comprising the exciton is shown in Fig. ||. 
The particles are confined approximately to one site for 
r > 0.4. For r < 0.4, they delocalizc rapidly, as indi- 
cated by the diverging radius. The self-trapped exciton 
evolves smoothly into a localized amplitude excitation of 



the charge-density wave, before disappearing into a band 
electron- hole pair at Tc{E), just as ^ the small polaron 
evolves through a large CDW-like polaron state before 
becoming a band hole at Tc{P)- 
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FIG. 5. Radius (r) of the exciton-hole and exciton-electron. 
A radius of corresponds to complete localization of the 
wave- function at one site. The radius diverges as F ^ Vc(E\ 
The finite value at Yc.^E) is caused by the finite system size 
of the numerical calculations. Inset: IPR-values of the gap 
states, indicating their degree of localization. They approach 
their asymptotic (symmetry determined) IPR values of i , i , 
\ (doubly degenerate), || and 1 for F ^ oo. 

IV. OPTICAL PROPERTIES 

In this section we discuss the observability of self- 
trapped excitons in experiment via an Urbach tail in op- 
tical absorption experiments and a multi-phonon spec- 
trum of the "breathing mode" in phonon experiments, 
both governed by the Franck-Condon effect. Zheng and 
Nasu have also obtained theoretically a tail in the op- 
tical conductivity. Their theory includes quantum lattice 
fluctuations (as does ours) and bears some relation to our 
theory. 

As a first step we neglect both Coulomb and lattice po- 
larization effects. Then the imaginary part of the dielec- 
tric constant can be computed exactly in electric dipole 
approximation. The dipole matrix element between an 
s state on a Bi atom at i? and one on the neighboring 
Bi atom at i? + ai is —imtax/Ti. This gives the correct 
intraband < k\p\k >= [m/h)de/dk relation between mo- 
mentum matrix element and group velocity. Then for for 
unpolarized light we get 

m'^uj-'Trri Jbz 
X 6(2^4 + -huj). (2) 

This was evaluated numerically using a tetrahedron code 
1^,^, and is shown as the dashed line in Fig.|[ The 
behavior of £2 is dominated by the divergence of the 
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joint-density of states at photon energy Ticu = 2A, cor- 
responding to transitions from states at the top of the 
lower Peierls band into those at the bottom of the upper 
band. These states are mainly responsible for the charge 
disproportion and have vanishing amplitude on the B 
and A sites respectively. These are also the states from 
which localized polaron and self-trapped exciton states 
form. Since transitions between these states are highest 
in optical density, we may expect the line shape to be 
dominated by the signature of self-trapped excitons. 

Materials with self-trapped excitons show Franck- 
Condon broadening of electronic transitions, a concept 
well known to molecular spectroscopists, since molecular 
excited electronic states generally have altered nuclear 
coordinates. The process involves simultaneous creation 
of electronic excitation and nuclear vibrations, causing 
the appearance of vibrational sidebands in the electronic 
spectra In simple solids, excitations are usually delo- 
calized which eliminates this effect However, in case 
of self-localized states, because of their significant lattice 
relaxation, Franck-Condon broadening of the electronic 
transitions must reappear. e2(w) in Franck-Condon ap- 
proximation is given by: 



e2iuj) oc \Px\' 



where Pw is the electronic dipole matrix element be- 
tween the electronic ground A and excited state A', w\n 
the probability to find n vibrational quanta in the elec- 
tronic ground state and {Xn\X'n') the overlap integral be- 
tween the vibrational wave functions XAn- This integral 
can be large even when the vibrational quantum numbers 
n, n' are quite different from each other, because the vi- 
brational wave- functions x\n, Xx'n' have origins which 
are off-set from each other. 

At zero temperature, the broadening of the electronic 
transition derives from the overlap of zero point motion 
of the oxygens with excited vibrations in the electronic 
excited state. The sum over n in Eq.(^ then has only 
one term, w\n = 5nfl . We lift the adiabatic approx- 
imation applied so far by adding quantized vibrations, 
but treat the problem in the atomic limit (t = 0), i.e. 
in single site approximation. The excited electron states 
couple to eleven oxygens, which are treated as indepen- 
dent Einstein oscillators with frequency loq and oxygen 
mass M in our model. They move in harmonic potentials 
and thus the Franck-Condon factor ||l2| is a product of 
eleven vibrational overlap integrals of the form |(0|n)p 
with: 



distribute n quanta over the eleven oxygens. The result 
is 
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(5) 



In our simplified Einstein model, absorption occurs at 
discrete frequencies corresponding to different numbers 
of vibrations created in the absorption process. In real- 
ity these peaks are broadened by phonon dispersion, and 
discreteness of the individual peaks should only be seen 
for one or two vibrational quanta. The line shape should 
look similar to the envelope function of Fig.^, with max- 
imal absorption at to — 2A and significant absorption 
below the gap at 2A=2 eV. 



w^n\{Xn\X'n')\H{niu + E^,, - Ey„,), (3) 
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FIG. 6. Imaginary part of dielectric constant versus fre- 
quency. The dashed line shows Eq.(^, where the integral 
was evaluated with a tetrahedron code [^,^. Taking lattice 
relaxation into account the absorption spectra will be domi- 
nated by the Franck-Condon effect Eq.(y), leading to signif- 
icant peak broadening and absorption in the gap, known as 
an Urbach tail. 

The absorption spectra of self-trapped excitons will 
change with temperature due to the temperature- 
dependent occupation of vibrational states in the elec- 
tronic ground state given by w\,n — e'xjp{l3nhojo) / Z , 
where the partition function Z is (1 — exp(/37ia;o))~^^ E^nd 
(3 = l/ksT . Temperature enhances the absorption in the 
low energy regime and shifts the on-set of absorption to 
smaller energies. The temperature effect up to room tem- 
perature is very weak. At T = 300K the probability to 
find a phonon on one of the eleven oxygens is less than 
6%. The analytical calculations for temperature correc- 
tions become increasingly tedious and Fig.^ shows a first 
temperature correction obtained numerically. 



a- \- [R- R'Y 



exp(--(i?-i?')'), (4) 



where a — y^Mujo/h and R and R' are the oxygen equi- 
librium positions of the ground and excited state re- 
spectively. The total absorption at a frequency hto — 
Ecxc + nhujQ^ requires a sum over all different ways to 
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FIG. 7. Imaginary part of dielectric constant (logarithmic 
scale) versus frequency at T=0, and the first temperature 
correction to the Franck-Condon spectra at 300K. 

The clearest evidence for self-trapped excitons is nor- 
mally 1^ observation of weak sub-gap features in opti- 
cal absorption a{oj) oc e2{u>) and a large Stokes shift in 
the luminescence spectrum. We are not aware of lumi- 
nescence measurements. Transient absorption measure- 
ments jl^ have been interpreted as indicating that inter- 
band excitations relax rapidly to a state shifted down by 
0.7 eV. Transmission measurements through thin single 
crystals were reported by Kozlov et al. Q. These have 
been interpreted as showing an indirect band gap at 0.5 
eV with anomalous temperature dependence. We believe 
that the measured absorption can be re-interpreted as an 
Urbach tail arising from the self-trapped exciton, with 
binding energy Coxc enhanced beyond our values by the 
on-site Coulomb U. 

Clear experimental support for the existence of self- 
trapped excitons in BaBiOs comes from harmonics of 
the Peierls breathing mode phonon seen both in Raman 
experiments [D and in the infrared These spectra 
are dominated by a strong peak at 570cm^^, assigned to 
the breathing mode Q — Aig phonon. Unlike in ordi- 
nary semi-conductors, a series of higher harmonics is ob- 
served. The simplest explanation for the Raman process 
is that the incoming photon creates a virtual self-trapped 
exciton, except that because of the Franck-Condon prin- 
ciple, the lattice coordinates stay at their ground state 
rather than their excited state values, which means that 
a cloud of virtual local vibrations is also created. When 
this virtual excitation re-radiates the scattered photon, 
not all the virtual vibrational quanta need to disappear. 
The test of this mechanism is that this process is res- 
onant when the laser photon energy coincides with the 
Peierls gap. This resonance was seen by Tajima et al. 

Similar behavior occurs in the Raman spectrum of 
LaMnOa ||T^ where a somewhat different species of self- 
trapped exciton is predicted to occur | [T5[ | . In the infrared 
spectrum of BaBiOa, the selection rule against Raman 
phonons appearing in infrared must be violated by lo- 
cal loss of inversion symmetry due to static or dynamic 
lattice fluctuations. Multiples of the breathing phonon 
are as easily excited as a single phonon because of the 
vibrational overlap factors. 
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